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We compare the equilibrium solution for the condensate 
obtained in the mean-field approximation to the master equa- 
tion for sympathetic cooling with the one obtained by Scully 
for a system in contact with a heat bath with the help of an 
analogy with the laser. While the mean-field approach yields 
analytical formulas for the approach towards equilibrium and 
for the equilibrium solution, it neglects the correlations be- 
tween occupation numbers of different single-particle states 
which are approximately kept in Scully's approach. Such ne- 
glect is admissible as long as the fraction of Bosons in the 
condensate does not exceed a few percent or so. 



Motivation. The experimental realization of Bose- 
Einstein condensation Q calls for a thorough theoretical 
understanding of the properties of the condensate. An 
important step in this direction was recently taken by 
Scully et al. [@||. Using the analogy between the con- 
densate and a laser, these authors derived the equilibrium 
form of the reduced density matrix po of the condensate 
for a system containing a fixed number Na of Bosons. 
The result is obtained as the steady-state solution of the 
time dependent non equilibrium master equation for an 
ideal Bose gas in a three-dimensional harmonic trap cou- 
pled to a thermal reservoir. The particle number Na 
serves as an important constraint. The solution found in 
Refs. [001 has the form 
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Here II" = \N)a o(N\ is the projector onto the lowest 
single-particle state containing N < Na Bosons. The 
temperature of the reservoir is denoted by T, while T c is 
a suitably defined transition temperature, and Zn is a 
normalization factor. 

Another starting point towards understanding the 
properties of the condensate is the master equation for 
sympathetic cooling. This master equation was, in very 
general form, derived in Ref. M. The cooling agent is a 
cooled gas in thermal equilibrium. The cooling mecha- 
nisms for a system in contact with a reservoir and for 
a system subject to sympathetic cooling W differ. In the 
first case, energy is exchanged via exciting or de-exciting 



the reservoir. In the second case, cooling is due to two- 
body collisions between the atoms in the cooling gas and 
those in the system of Na Bosons. The form given in 
Eq. for the density matrix of the condensate is also 
obtained from the master equation for sympathetic cool- 
ing (|] by postulating that the elements of the density 
matrix referring to excited single-particle states have at- 
tained the equilibrium form but that this form depends 
on the number N of Bosons in the ground-state con- 
figuration \N)q. The agreement between the forms of 
Po found from two different approaches gives confidence 
that Eq. ([l]) represents a good approximation to the true 
from of the density matrix for the condensate. 

Recently, the present authors have developed a mean- 
field approach to the master equation for sympathetic 
cooling Q . This approach yields explicit analytic expres- 
sions for the density matrix of both, the condensate and 
the A^-Boson system in excited single-particle states. 
The approach yields not only the equilibrium solution 
but also the time-dependence of the cooling process. The 
reason is that the mean-field equations possess a SU(1, 1) 
dynamical symmetry. Therefore, the equations are inte- 
grate 0j§|. 

In the present Letter, we compare the equilibrium so- 
lution for the density matrix pa of the condensate as ob- 
tained from the mean-field approximation with Eq. (Q). 
The comparison will cast new light on both, the form 
of po and the limitation of the mean-field approach. In 
addition, we briefly display the time evolution towards 
equilibrium as obtained in the mean-field approach. 

Master Equation. Starting point is the master equa- 
tion derived in Ref. . We use the notation of Refs. j|,U ■ 
We label the system subject to sympathetic cooling as 
system A. This system consists of Na Bosons. The mas- 
ter equation for the dependence of the reduced density 
matrix pA{t) for system A on time t reads 



dp A (t) 
dt 
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Here, Ha is the sum of the single-particle Hamiltonians 
for the atoms in system A (each containing the kinetic en- 
ergy and the harmonic trap potential) while H' A _ A repre- 
sents the (weak) interaction between the atoms in system 
A. This interaction is neglected in what follows because 
sympathetic cooling is used precisely when H' A _ A is very 
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small. The action of the Liouvillean C on the reduced 
density matrix pA{t) is given by 



n : n' ,rh,rh' 
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The single-particle states of the three-dimensional 
isotropic harmonic trap are labelled by the quantum 
numbers rh = (m x ,m y ,m z ) with m x ,m y ,m z integer. 
The creation and annihilation operators for these states 
are written as at^ and a^a, respectively. The index is 
used for the non-degenerate ground state of the trap. 
The rate coefficients r™^ 1 are given in Ref. 0. We do 
not repeat the definition here. Suffice it to say that these 
coefficients account for the interaction between particles 
in system A and those in the cooling system B. We as- 
sume that the number Nb of atoms in system B is very 
large. Then decoherence acts very quickly |^] and reduces 
the density matrix to diagonal form. Hence, 

cpA = E r s( 2a ii a «p.4W4 a ™ 

-ata^atanPAit) - PA(t)aia A a'laf i ) . (4) 



Mean-Field Approximation. Eq. (jj) serves as the 
starting point for the mean-field approximation. We ap- 
ply this approximation in standard fashion by replacing 
on the right-hand side of Eq. (||) one pair of creation 
and annihilation operators referring to the same single- 
particle state rh or h by its expectation value (a^a^). 
The quantity 
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is the mean particle occupation number of the single- 
particle state rh. For the term CpA, this procedure yields 

CpA = J2 T t'ri. N * l 2a lnPA(t)arn 

rh^n 

-a^a^pAit) - pAitfa^at - p A ^j 
+ E r ^ Wn + V (^PA(t)al n 

rh 

-a^arnpAit) - pA^a^a^ + p A ^j ■ (6) 

Central to the mean-field approach is the assumption 
that there are no correlations between the occupation 
probabilities of the harmonic trap levels n and rh. 

Rate Equation. Eqs. (|J) and (||) can be reduced to a 
rate equation. We use a basis in Hilbert space defined 



by a product of all single-particle states m, each such 
state being occupied by Bosons. We take the trace 
of pA{t) over all single-particle states n with h ^ rh 
(this includes a summation over all occupation numbers 
Nfij) and denote the result by p^i(i). That same notation 
was used already in Eq. ([[]) with rh = 0. We recall that 
PA{t) is diagonal in energy representation. It follows that 
Pm(t) can be written as a sum of the projectors 11^ = 
\N) AA (N\, 



N A 

Pra(t) = E P > 
JV=0 



N 



(*) nil 



(7) 



The time dependent mean field occupation probabilities 
P£(t) differ from zero only for N < Na- Taking the 
corresponding trace of Eq. (||) and using Eq. (||), we find 
that the P^(t) obey the rate equation 
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The cooling and heating coefficients K m and H m are 
given by 
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^From Eqs. (^|,|7|) and (g) we obtain for the mean occu- 
pation number iV,%(i) 
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This equation is consistent with conservation of particle 
number, J^a, d-N^/df = 0. 

Condensate. Putting rh — in Eq. (|8|), one obtains the 
equation for the probability distribution of the number 
of atoms in the single-particle ground state, i.e., the con- 
densate. The equilibrium solution obeys dP^(t)/dt = 0. 
It is easily seen that the resulting equations for the time- 
independent coefficients P^ do not possess a non-trivial 
solution. This is because the constraint = for 
N > Na is too rigid for the mean-field approach. We 
relax this condition, write 
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and require only that the sum ^2^ =Na P^ be negligibly 
small. The resulting equation for the P^'s reads 

Vj 
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The solution is 



plV _ pO 



'o where x — 



(13) 



Eq. (|T^) implies that in the stationary case we must have 
X < 1. The normalization condition yields Pq = (1 — x). 
The mean value Nq is given by Nq — x/(l — x)- Con- 
versely, we may replace x everywhere by iVo / (Nq + 1 ) . To 
discuss the validity of the mean-field solution, we impose 
the constraint that En=n a P o l P o = X Na /(1 - x) < 
exp(-a). For N A > 1, this yields N /N A < 1/a. This 
condition applies as long as a > ln(iVo) and shows that 
for the condensate, the mean-field approximation (as de- 
fined in the framework of this paper) begins to fail when- 
ever the ratio N /N A grows beyond a few percent or so. 
Similar conclusions apply to the stationary occupation 
probabilities Pj~ of excited single-particle states, except 
that here we do not expect the mean occupation number 
ever to approach values close to Na- 

The Pq' ~ x N given by Eq. ( |l3| ) decrease monotoni- 
cally with increasing N. This behavior is in marked con- 
trast to that of the equilibrium solution Eq. p) found 
first by Scully and reproduced, in the present con- 
text, in Ref. jq]. The ground-state distribution as given 
by Eq. (|l|) is Poissonian, except for the additional nor- 
malization factor. 

The physical situations encapsulated in Eqs. (Q) and 
( |l3| ) are completely analogous to the behavior of a laser 
far above and below threshold, respectively M. The dif- 
ference between Eqs. (Q) and ( |l3| ) is due to the strong 
correlations between the numbers of particles in excited 
single-particle states and in the ground state. That 
correlation is taken into account approximately in the 
derivation of Eq. (fil). Indeed, in this derivation it is 



assumed that for m / we have (a^a^) 



with Nfji(N) constrained by N, the number of Bosons 
in the ground-state configuration \N)q. The resulting 
nonlincarity is analogous to the nonlincarity which gov- 
erns the behavior of a laser far above threshold. It carries 
the solution beyond the mean-field approximation and is 
well suited to describe the probability distribution of the 
atoms in the condensate whenever their number is of the 
order of Na- As mentioned in the beginning, there is 
good reason to believe that Eq. ([I]) correctly describes 
this situation. In the mean-field approach, on the other 
hand, the correlations between the numbers of particles 
in the excited single-particle states and in the ground 
state are totally neglected. This situation is analogous 
to the behavior of a laser below threshold (linear case). 
We conclude that Eqs. (|l|) and (|l3|) describe two differ- 
ent regimes. The equilibrium solution Eq. ( |l3| ) describes 
a system where at most a few percent of the atoms are 
in the ground state. As that fraction increases, there 
is a gradual transition to another regime described by 
Eq. (Q). The intermittent case is apparently not covered 
by either formula. 



In the remainder of the paper, we demonstrate that 
the mean-field approximation is well suited to describe 
the approach of the system towards equilibrium, and not 
only the stationary case. This is due to the fact that 
the master equation is linear in the generators of the 
group SU(1, 1). Details of the calculation are given else- 
where Writing the time-dependent master equation 
in the form 
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with r(t) dependent upon time, we construct a time- 
dependent similarity transformation T(t) which diago- 
nalizes T(t). Because of the form of T, there exist two dif- 
ferent similarity transformations which accomplish this 
aim. However, only one of the two fulfills the condition 
for the viability of the mean-field solution, namely, that 
the coefficients multiplying 11^ vanish asymptotically for 
large N. Using this transformation, the time-evolution 
of the reduced density matrix can be determined and is 
given by 

p A (t) = exp[a+(t)^+] exp[a-(t)K~] 
xexpQ^ dr[j(T)K - K m (r) - H m (r))] \ PrA (0) (15) 

where Pm(0) is fixed by the initial condition, and j(t) = 
AH m {t)a+{t) - 2[K m (t) + H m (t)]. The time-dependent 
functions a± (t) are the solutions of the differential equa- 
tions 
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The initial conditions are a±(0) = 0. The action of the 
operators K + , K~ and K° on the quantities n^ fc = 
\n)rnrh(k\ is given by 
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^From the asymptotic behavior of the coefficients K m , 
H m , 7 and a±, it can be shown that as t — > oo, Pm(t) 
approaches the equilibrium solution which for fh = is 
displayed in Eq. (|l3|). 

In summary, our results show that the mean-field ap- 
proach is a useful tool to investigate both sympathetic 
cooling and properties of the condensate. The limiting 
factor in the approach is the neglect of correlations be- 
tween occupation numbers in the ground state and in 
excited single-particle states. Such correlations generate 
a non-linear term in the cooling and heating coefficients 
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and become important whenever the fraction of Bosons 
in the condensate exceeds a few percent or so. While the 
master equation for the mean-field approach allows for 
an analytical solution, it is unlikely that such a solution 
exists in the general case. The distributions of occupa- 
tion probabilities of states with No Bosons in the ground 
state differ very much in the mean-field approach and 
in the approximation obtained by Scully. The latter is 
probably adequate whenever the fraction of Bosons in 
the ground state is of order unity. The two distributions 
display similarities with the behavior of a laser below and 
far above threshold, respectively. 
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